Experimental Investigation of Plastic Deformations Before Granular Avalanche 
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We present an experimental study of the deformation inside a granular material that is progres- 
sively tilted. We investigate the deformation before the avalanche with a spatially resolved Diffusive 
Wave Spectroscopy setup. At the beginning of the inclination process, we first observe localized 
and isolated events in the bulk, with a density which decreases with the depth. As the angle of 
inclination increases, series of micro- failures occur periodically in the bulk, and finally a granular 
avalanche takes place. The micro-failures are observed only when the tilt angles are larger than a 
threshold angle much smaller than the granular avalanche angle. We have characterized the density 
of reorganizations and the localization of micro-failures. We have also explored the effect of the 
nature of the grains, the relative humidity conditions and the packing fraction of the sample. We 
discuss those observations in the framework of the plasticity of granular matter. Micro- failures may 
then be viewed as the result of the accumulation of numerous plastic events. 

PACS numbers: 83.80.Fg,47.57.Gc,45.70.Ht 



I. INTRODUCTION 

The onset of avalanches in granular materials is a well- 
documented area of research motivated by evident appli- 
cations in the field of risk prevention. A typical granular 
experiment to study such phenomenon consists in the 
slow inclination of a sand filled box until the destabiliza- 
tion of the sand pile p]-[3| . A lot of studies have been de- 
voted to the angles at which this avalanche initiates and 
stops @,[5|. The critical angle, or avalanche angle, is then 
defined as the angle at which a large flow of grains oc- 
curs, with the surface of the heap stabilizing at a smaller 
angle, the angle of repose. Usually studies of stability of 
granular heaps focus on the succession of avalanches 
either in rotating drums or by using a local perturbation 
of a heap. In those two cases, the heap has been built by 
successive avalanches, giving a specific internal structure 
to the bulk. The time intervals between successive events 
and the size of the avalanches is then discussed, one of 
the points being to know if those distributions are or 
are not power-law distributions as predicted for a critical 
phenomena ^ . In the case of a granular pile without ini- 
tial inner structure progressively inclined in the gravita- 
tional field, specific types of rearrangements taking place 
before the first avalanche have been identified. Experi- 
mentally, two kinds of movement in the heap previous to 
its failure have been detected: small rearrangements, im- 
plying only a fewgrains, and periodic precursors of large 
amplitude 0-EJ- Most of the studies reporting this 
phenomenology are based on the observation of the free 
surface of the pile, but recent works give clues that these 
precursors are bulk phenomena 0, 0. To our knowledge, 
these regular events have never been reproduced numer- 
ically or explained theoretically, although preavalanche 
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instabilities showing an intermittent evolution have been 
studied numerically [To|, [ll| . 

The periodicity of the precursors are reminiscent of the 
stick-slip phenomenology, which is usually invoked as an 
explanation of the observed regularity. To our knowledge, 
a complete interpretation in that framework, including 
the prediction of the starting angle of the process, which 
is well beyond the angle of friction of the material, or 
of the selected period, has never been made. Stick-slip 
behavior emerges when a frictional system is submitted 
to an increasing load at a small enough rate [12| . It orig- 
inates from the difference between the static friction and 
the dynamic friction. In the simplest model, a traction is 
exerted on a frictional slider by a spring. Elastic energy 
is first stored in the spring while the slider is stuck until 
the tangential force is large enough to overcome the static 
friction. The slider then enters an oscillating response 
that leads it to stick again when the velocity vanishes. A 
model inspired by this principle has been used to explain 
the succession of avalanches in a rotating drum [12j. In 
the same spirit, arrays of sliders connected by springs 
(Burridge-KnopofT model) [l3| are used to model seismic 
processes. Fault gouges are often modeled in the geology 
literature by a granular material sheared between two 
planes and the stick-slip behavior of such systems has 
been extensively studied in an established regular regime 
of shearing of the granular layer [14], ll5| . The existence 
of periodicity in seismology is still under discussion and 
recent models exhibit the coexistence of quasiperiodicity 
and criticality [l6| . 

It is striking to observe that a very similar phenomenol- 
ogy of regular precursors of rupture before the establish- 
ment of the regular stick-slip response can be found in 
the literature in two other configurations closely linked 
to the problem of failure in granular materials: the on- 
set of frictional motion and the mechanical response of 
metallic glasses. In the case of frictional motion, regular 
precursors are observed before the onset of sliding 
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[20| . The present interpretation relies on the inhomo- 
geneity of the spatial repartition of the normal and shear 
stresses at the interface leading to the possibility of reach- 
ing the Coulomb criteria in some localized areas before 
that the whole system is destabilized [U, [22| . The links 
with the models of earthquakes described in the previ- 
ous paragraph have been made [23|, [24|. In the second 
system, the case of the mechanical response of metallic 
glasses, the observed stick- slip phenomenology is called 
serration. The mechanical response of such amorphous 
glasses is very close to the response of granular mate- 
rial and displays similar features: creep, shear-banding 
and precursors to the rupture, all of which seem to be 
generic to amorphous materials. Comparing the loading 
curves obtained for metallic glasses (Figure 1(a) in |25j), 
granular materials (Figure 2(b) in [26]) and onset of fric- 
tion (Figure 1(a) in [l7[) a striking similarity is observed. 
Each of these three curves displays small regular stress 
drops during the otherwise quasi-elastic loading until the 
failure (or sliding) of the material. Those drops begin 
to appear well before the rupture of the material or the 
onset of motion of the slider. After the yield, the system 
enters a regular stick-slip motion. In the case of gran- 
ular materials, it has been shown that the value of the 
stress at the first micro-rupture event is linked to the 
stress at which the global yield will take place [26| : inde- 
pendent of the packing fraction of the system, the ratio 
between those two stresses is constant. Using a method 
of detection identical to the one that will be used in this 
paper, it has been shown that the drops in the loading 
curve are linked to the appearance of a shear band in the 
system [27j. During the loading of the system, "spots" 
of localized strong deformation are also observed. As it 
will be shown, these localized spots tends to appear at a 
higher rate, and to cluster, just before a shear band event 
occurs [27^ 

Consequently, a natural question emerges regarding 
the understanding of the nature of those precursor events 
and their potential link to the small rearrangements. 
With this aim in mind, in this article, we study ex- 
perimentally the response of granular material to a pro- 
gressive, quasi-static, inclination. We will focus on the 
description of the behavior of the heap before the first 
avalanche. For the study of the rearrangements and of 
the precursors, we use an original method of measure- 
ment of small deformation based on Diffusive Wave Spec- 
troscopy [28|, [29| . This method gives access to the phe- 
nomenology of the rearrangements and the precursors at 
the side of the sample, and consequently throughout the 
depth of the pile, giving hints of what happens in the 
bulk of the system. 

The structure of the article is as follows: in part ^TT] 
we will describe the mechanical setup ( £111 Aft , the in- 
terferometric method of detection of the rearrangements 
and deformations ( £111 Bp , and the granular materials and 
the experimental protocols of preparation of the sam- 
ples ( £111 Cp . In the third section we will present the ob- 
served experimental behavior. After an overview of the 



experimental observations ( £1111 Aft , we will characterize 
the periodic precursor events ( £1111 Bp and then the small 
rearrangements phenomenology ( £1111 Cp . We also discuss 
the influence of several parameters (shape of the grains, 
humidity and packing fraction of the sample) on the ob- 
served phenomenology ( ^III Dp . The results are discussed 
in £1IV1 In £1IV Al we compare the observations with the 
classical Mohr- Coulomb theory of failure. In £1IVB[ we 
interpret these events in the framework of plasticity of 
granular materials, and coupling between plastic events. 



II. EXPERIMENTAL SET-UP 

A. Mechanical setup 

The mechanical set-up consists of an optical board that 
can be inclined continuously by a crankshaft system (see 
Fig.[]Ja)). The laser source that illuminates the sample, 
the optical detection setup and the granular container 
are all fixed on the board so that the whole set-up rotates 
with the granular pile (see Fig.[T]). The translation of the 
crankshaft is carried out by a motorised linear translation 
stage (Micro Controle UT100-50PP) controlled electroni- 
cally (Micro Controle ITL09). This system ensures a low 
level of vibration at a small rotation rate. For transla- 
tions between 25 /xm/s and 100 /xm/s, the tilt rates range 
between 0.02° /s and 0.08° /s. We will show in the follow- 
ing that all the experiments can be considered to be done 
in a quasi- static limit. 




FIG. 1. (a) Schematic of the crankshaft system with the tilt- 
ing optical board, (b) Schematic of the setup clamped to the 
tilting optical table. The laser source is expanded in order 
to illuminate the side of the sample. The backscattered light 
is collected on a CCD camera. A lens images the side of the 
sample on the CCD and a diaphragm controls the size of the 
speckles on the camera. 

The granular sample is contained in a rectangular box 
of dimensions 20 cm x 14 cm x 8 cm, the lateral sides of 
which are glass windows. One lateral side is illuminated 
by a laser (Melles Griot, 15 mW, 633 nm). The beam 
of the laser is expanded by a lens in order to illuminate 
the entire surface of one side of the box. A lens images 
that side of the sample on a CCD camera (Prosilica GC 
2450, 2448x2050 pixels, cell size 3.45x3.45 /mi), with 
a magnification 7 = 0.09. Because of the disorder of 
the granular media, the backscattered light collected has 
a speckle pattern in intensity. A diaphragm allows us 
to adjust the size of these speckles on the CCD camera 
(see next section for the optical method). The camera 
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is interfaced with a computer and a Labview program 
ensures the acquisition of successive images at a frame 
rate that has been chosen between 1 fps and 3 fps in the 
experiments shown here. 

B. Optical method 

The measurement of deformation and rearrangements 
in the granular material is carried out using an in- 
terferometric technique based on Diffusive Wave Spec- 
troscopy (DWS). The method has been extensively de- 
scribed in [28|, [29| . The treatment of raw experimental 
images is based on the correlation of backscattered in- 
tensity between two successive images following a multi- 
speckle procedure. The link with the corresponding de- 
formations that have occurred inside the sample is made 
using a model [lH, [30| . 

The schematic of the optical part of the setup is shown 
in Figure QJb). The beam of a laser source is expanded 
to illuminate the granular sample and the side of the 
sample is imaged on a camera with a lens. Because of the 
disordered structure of the scattering material, light rays 
inside the media follow random paths. The backscattered 
rays explore a typical zone of length I* in the sample, 
called the transport mean free path of light. For glass 
spheres of diameter d, I* ~ 3.3d [30|. As the incident 
light is coherent, those backscattered rays interfere and 
give a speckled structure to the collected intensity. The 
size of the coherence areas of the speckle in the image 
plane is controlled with a diaphragm. When the sample is 
deformed, the scatterers are displaced and the light paths 
modified so that the interference pattern changes. These 
changes can be quantified by calculating the correlation 
of intensities between two states. 

In a backscattering configuration, a light ray typically 
explores a volume of characteristic size (P) 3 inside the 
scattering material. By using domains of size (7/*) 2 in 
the images for the ensemble average, we can then obtain 
a local evaluation of the correlation function with the op- 
timal spatial resolution. Correlation maps are thus ob- 
tained, giving insights into the processes that have taken 
place in the sample during the loading. 

The intensity correlation between two successive im- 
ages, labelled 1 and 2, is calculated using the relation: 

{hh) - (h)(h) 

1 V(A 2 )-(A)V(^)-(^> 2 ' 

where I\ and I2 are the scattered intensities on images 
1 and 2. The electronic noise of the camera and possi- 
ble stray light are taken into account with a procedure 
explained in [31]. The ensemble averages (•) are calcu- 
lated over square sub-areas of the images of size 16x16 
pixels, which typically contain about forty coherence ar- 
eas. From initial images of size 2448x2050 pixels we then 
obtain maps of 153x128 pixels. With the lens magnifica- 
tion chosen here, the spatial resolution of the final map 
is about 608 fim. 



The map of correlations are calculated between suc- 
cessive images, giving information about the incremen- 
tal decorrelations that have taken place in the system 
for an angle increment. When the tilt rate was modi- 
fied, the frame rate of the camera was chosen to ensure 
that the angle increment between two successive images 
stayed between 0.06°/slice and 0.08°/slice. This allowed 
for consistency and ease of comparison between the dif- 
ferent experiments. 

In the case of linear elastic samples [28], the deforma- 
tion of the material corresponding to the measured decor- 
relation can be estimated using a model. The extension 
of that model in the case of granular materials has been 
discussed in [29|. Considering that the deformation be- 
tween the two correlated images is affine at the scale of I* 
in the material, that the light is strongly scattered in the 
sample, and that we observe the backscattered light, the 
expected dependence of the intensity correlation function 
on the strain tensor is as follows: 



Gj-expt-cvTM), (1) 

where f (XJ) is a quadratic function of the strain tensor U, 
(see [28|,l29| for details of that function), corresponding to 
the deformation increment that has taken place between 
the two states associated with the images used for the 
calculation of the correlation [28|, [29|. Following [29j . 

the constant c has been taken as c = 2?7 2 ^ L /*^/|, with 
77 ~ 2 a constant depending of the optical setup and A 
the wavelength of the laser. The typical deformations 
probed by this method for values of the transport length 
similar to the one of the glass beads used here, I* ~ 
990 fim, are between 10~ 6 and 10~ 4 HzlH^. For sand, 
we measured I* ~ 660 jim c± 2d. This lower value of l*/d 
compared to glass beads may be understood by noticing 
that the surface of sand grains should itself diffuse light. 
However, since the ratio l*/d is not very different for 
glass beads and sand, we expect the magnitude of probed 
deformations to be very similar. 



C. Granular samples and preparation 

Experiments have been performed for two different 
type of grains: glass beads (Silibeads, Figure Ufa)) and 
sieved sand (Sifraco, Figure EJb)) with a similar range of 
diameters: 200-400 /am as can been seen on the cumula- 
tive size distributions for the two materials in Figure [2j 
We observe in the insert pictures in Figure [2] that the 
glass beads have a lot of defects and that the main dif- 
ference between the two kind of grains is that the sand is 
faceted whereas the glass beads are rather spherical. The 
angle of avalanche for the glass beads is about 29° and 
for sand about 35.5°. The cumulative distributions show 
the presence of few particles in the range 50-100 jim for 
the two types of granular material. Since the dynamics of 
precursors is in agreement withpreviously reported ex- 
periments with larger beads 0, [?H2| , we may think that 



the details of the particle size distribution is not crucial. 










100 200 300 400 500 600 700 800 100 200 300 400 500 600 700 800 

beads diameter (jim) beads diameter (|jm) 

FIG. 2. Cumulative distribution of the grain diameter for the 
glass beads (a) and the sand (b). Pictures of the granular 
materials are shown in inserts. The median size of the grains 
is about 280 /xm for the glass beads and 330 /im for the sand. 

Three different levels of compaction have been used. 
For the lowest compaction samples, a grid (mesh size 
4 mm) is placed at the bottom of the empty box, the 
granular material is gently poured in the box, and then 
the grid is slowly lifted through the sample. More com- 
pacted configurations are obtained from this initial loose 
preparation by tapping the filled box in order to increase 
the density of the system. Grids of different mesh sizes 
and forms have been tested not to change the global phe- 
nomenology. The less compacted samples have packing 
fraction around C — .58 ± .02, the intermediate case 
(called "normal" in the following) has C ~ .61 ± .02, 
and the highest packing fraction is close to G ~ .64. 

Four different relative humidity conditions have been 
used. Ambient humidity was between 25 % and 45 % 
(called "normal" condition in the following). We obtain a 
humidity smaller than 20 % by leaving the granular sam- 
ple in an oven overnight and then letting it cool down 
in presence of desiccant. We obtain high humidity by 
leaving the granular sample overnight in an enclosure in 
the presence of either saturated salt water or pure water. 
The relative humidity in the enclosure is controlled by 
the nature of the salt, and we then obtain samples pre- 
pared at respectively 70 % and above 80 % humidity. For 
all those values of the relative humidity, the dynamics is 
unchanged when modifying the duration of the experi- 
ment. So we can neglect relative humidity variations at 
the timescale of the experiment. 

III. EXPERIMENTAL RESULTS 

A. Main phenomenology 

A typical movie of the observed behavior in the so- 
called "normal" conditions (packing fraction around .61 
and relative humidity of 25 %) is given in Supplemental 
Material [32j. The granular material used in this movie 
is glass beads but the phenomenology is very similar to 
sand. Each image is a map of the correlation between 
two successive images and corresponds to the incremen- 
tal deformation during an angle variation of 0.08°. The 
color scale is the following: white correspond to Gj = 1 
(maximal correlation) and black to Gj = (vanishing 



correlation). The value of the angle of inclination of the 
board is indicated in degrees at the bottom-left of the 
film. The free surface of the granular pile can clearly 
be identified due to the fact that the light that does not 
come from the sample is totally decorrelated. The size of 
the area seen in the film is 7 cm x 8.5 cm from the upper 
middle of the box (see the schematic in Figure [3]). In or- 
der to evidence the small rearrangements, the contrast of 
the maps shown on Figure [3] has been enhanced using a 
threshold value of Gj = 0.75, so that all the correlations 
smaller than that value appear in black. 



free surface — > 




8.5 cm 



FIG. 3. Schematic of the area used in the movie in Sup- 
plemental Material [33]. Two examples of correlation maps 
extracted from the film, for angles 10° and 27°, are shown 
and exhibit typical rearrangements. The decorrelation maps 
of these images have been thresholded for Gi smaller that 
0.75 to allow the rearrangements to be clearly identified. The 
granular material is glass beads, the packing fraction around 
.61 and hygrometry of 25 %. The convention used for the 
axes are indicated on the figure. 

Two distinct types of phenomenology can be observed. 
Firstly, spots of different sizes appear well before any 
macroscopic event in the sample. Such events are shown 
for example on Figure [3] for the 10° angle. Enhance- 
ment of the contrast shows that such events appear as 
soon as the inclination process begins. These spots are 
seemingly of the same nature as the rearrangements de- 
scribed previously in the literature [3| and which were 
observed at the free surface of the sample. Our observa- 
tion tends to show that those rearrangements happen in 
the whole bulk and are not mere displacements of grains 
due to irregularities of the free surface originating from 
the preparation of the sample, as has been postulated in 
previous studies. The density and the intensity of those 
events increase with the angle. The density seems also to 
depend on the depth in the sample. Such localised events 
are reminiscent of the spots observed in a shear granu- 
lar sample in [27]. The fact that such spots have been 
observed at different borders in several systems limited 
by different boundary conditions supports the hypothe- 
sis that such localized events happen in the bulk of the 
sample. A quantitative estimation of the energy released 
locally compared to the overall work done during creep 
processes J27[ also supports that affirmation. The spots 
were identified in [271 as the localized events of deforma- 
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tion introduced by Argon [33j to describe the plasticity 
of amorphous materials. 

Secondly, from an angle of about 15°, large and almost 
regular events begin to happen, which appear as succes- 
sive large decorrelation zones. These events correspond 
to the periodic precursors already described in the lit- 
erature [3[. Our experiment makes it evident that these 
precursors involve a large part of the bulk. These pre- 
cursors begin to happen from about 15° and occur nearly 
regularly, with an angular periodicity between 1° and 2°. 
Quantitatively, these values correspond well to the values 
that have been reported in the literature for other sizes 
of grains and containers d, 0-0 • We can observe that a 
precursor mobilizes a slice of the sample parallel to the 
free surface. Such precursors are in fact micro-rupture 
events in the bulk material. We also observe that the 
position of the micro-rupture occurs at a depth that in- 
creases with each consecutive event. Between two pre- 
cursors, local rearrangements are also still observed. Be- 
cause of the approximative translational invariance of the 
phenomenon along the optical plane, we obtain a spatio- 
temporal representation by averaging the correlations at 
constant depth in each image along horizontal lines (x- 
axis): (Gi(0; + AO; x;z)) x . This gives a ID representa- 
tion of the dependence of the correlation with depth for 
each image. By juxtaposing these ID lines for successive 
angles of inclination, we obtain spatio-temporal graphs 
such as the one shown on Figure HJ As the dependence 
between time and angle is linear in a good approxima- 
tion, the horizontal axis can be considered as represent- 
ing either time or angle. On such a representation the 
periodic precursors are displayed distinctively while the 
rearrangements contribute to an average decrease of the 
correlation, giving a blurred or shaded aspect to the im- 
age where the activity in terms of small rearrangements 
is high. 

avalanche depth 
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FIG. 4. Spatiotemporal diagram obtained from the correla- 
tions map by spatially averaging each horizontal line of the 
images: (Gi(0; + AO; x; z)) x . The averaging procedure gives 
the mean decorrelation as a function of depth for each image. 
The granular material is glass beads, the packing fraction 
around .61 and relative humidity of 25 %. Correlations are 
calculated between images acquired at angles and + AO, 
with AO = 0.08°. 

It is noticeable in the movie or in Figure [3] at the angle 
of 27° that before and after a micro-rupture event oc- 
curs, a large amount of activity can be observed at the 



depth where the failure will take place. Larger spots of 
decorrelation are observed in the zone that will fail com- 
pared to those that appear at a larger depth. The decor- 
relation events seem to align and cluster at the place 
where the rupture will take place, reminiscent of the 
phenomenolo gy already observed in a sheared granular 



material [26|, [27j where regular micro-ruptures were also 



observed before the final yield stress of the material is 
attained. 

In the following section we give a more detailed study 
of that global phenomenology which is very robust and 
resistant to changes of parameters (grain shape, relative 
humidity, packing fraction). Firstly we will detail the 
precursors events, then we will describe the rearrange- 
ments. Finally, we will study the changes of the phe- 
nomenology with hygrometry and packing fraction. 



B. Precursors 



1. Periodicity 



We have studied experimentally the dependence of the 
periodicity with respect to the tilt velocity. Figure [5] 
shows the average period of the precursors in degrees for 
different rotation rates for the same conditions of relative 
humidity and packing fraction (so-called "normal" condi- 
tions: packing fraction about 0.61 and relative humidity 
about 30 %). Each point corresponds to the average pe- 
riod between the precursors during one run. The value of 
the period slightly decreases when the tilt rate increases, 
but overall the phenomenon seems not to depend on the 
tilt rate. This mean that the experiments can be consid- 
ered as having been carried out in the quasi-static limit. 
We also see in Figure [5] that the nature of the material, 
glass beads (green (light gray) points) or sand (black tri- 
angles), does not seem to modify significantly the overall 
mean value of the period which is a little over 1° for both 
materials. The error bars give the minimal and maximal 
values of the period during each run. That dispersion 
around the mean value follows a trend during one ex- 
periment: a significative increase of the period, given by 
the size of the error bars, has been observed between 
the first occurrence of the precursors around 15° and the 
avalanche angle. 

To conclude, the periodicity does not depend of the tilt 
rate in the quasi-static limit for which the experiments 
have been done. It is therefore justifiable to express that 
periodicity in terms of angles and not in terms of time 
intervals. The period has been found to be about 1° for 
the two types of granular materials. It has to be noted 
that similar values of the period have been obtained for 
much larger beads (2 to 3 mm), in containers of different 
sizes and materials d, 0-0 • 
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FIG. 5. (Color online) Average period of the precursors with 
the tilt rate of the sample. The error bars give the extremal 
values of the period around the mean value during one ex- 
periment. Green (light gray) points are for glass beads, black 
triangles for sand. The horizontal lines give the average value 
of the period over all the runs for each type of material. All 
the experiments have been done with a similar preparation: 
packing fraction about 0.61 and relative humidity about 30 % 



2. Size of the precursors 

An important question that remained unanswered in 
previous works was whether the bulk is involved in the 
rearrangements and precursors phenomena. As a mat- 
ter of fact, most of the previous experiments have been 
carried out using only free surface observations [H, Q , so 
that rearrangements and precursors could be supposed 
to be a mere free surface phenomenon. On the other 
hand, acoustic measurements 0, Hj gave rather indirect 
indications of the bulk mobilization. In our experiment 
we are able to visualize directly how deep the phenomena 
go inside the sample, even if the observation is confined 
to a thin layer near the wall of the container. We show 
that our results reinforce the hypothesis that a part of 
the bulk is indeed mobilized in a precursor event. The 
volume that is involved increases as the system comes 
closer to the avalanche angle. 

The successive depths of precursors can be obtained 
from the spatiotemporal graphs such as the one of Fig- 
ure [31 The values of the depth of the peaks, normalized 
by the diameter of the grains, as a function of the angle 
at which they appear are shown on Figure [6j The main 
part of the graph shows the evolution of those depths for 
different typical experiments. We can observe that the 
dependence of the size of the precursors with the inclina- 
tion angle is roughly linear, even if in some experiments 
brutal "jumps" can be observed (results displayed with 
solid lines). The green (light gray) lines correspond to 
experiments with glass beads and the black lines to the 
ones with sand. The slope of increase of depth with the 
angle is generally larger for the glass beads than for the 
sand. Also, the first precursor tends to appear sooner for 
glass beads than for sand. The insert of Figure [6] shows 
all the results from 25 experiments for glass beads and 14 
for sand superimposed, giving an idea of the dispersion of 
the results from one experiment to the other. All the ex- 
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FIG. 6. (Color online) Depth z, in units of grain diameter 
d, of the successive precursors plotted as a function of the 
angle of inclination of the granular sample 6. All the ex- 
periments have been done for "normal" compaction (~ .61) 
and relative humidity conditions (~ 30 — 40 %) but different 
tilt rates. Main graph: examples of evolution for glass beads 
(green (light gray) lines) and sand (black lines) showing the 
roughly linear dependence of the depth with the angle. In- 
sert: gathered results for a total of 25 experiments for glass 
beads (green (light gray) points) and 14 experiments for sand 
(black triangles). Lines are linear fits of those data: for glass 
beads, z/d = 9 x (0 - 14°), for sand: z/d = 4 x (6 - 16°). 

periments have been done in the "normal" conditions of 
relative humidity and packing fraction and for tilt rates 
between 0.02° /s and 0.08° /s. Linear fits using the whole 
set of data for each type of granular material have been 
done and are shown as solid lines in the insert. The slope 
of the line is about 9 beads diameter per degree for the 
glass beads and 4 for the sand. For the two types of ma- 
terial, the intersection of the linear fit with the abscissa 
axis is about 15°, even if the first precursors appear later 
in sand than in glass beads. 

Thus, precursor events involve a fraction of the ma- 
terial that increases roughly linearly with the angle of 
inclination. That increase is approximatively two times 
larger for the glass beads than for the sand. The part 
of the sample that is mobilized is parallel to the surface, 
such that the process seems to imply a rectangular block 
of granular material in the upper part of the sample. A 
closer observation of the movie in Supplemental Mate- 
rial [32| shows that the surface of separation between the 
mobilized zone and the quiescent zone bends toward the 
bottom of the box, which could be an effect of the side 
boundaries. A more detailed description of the processes 
happening in the decorrelated area during a precursor 
event is difficult giving the impossibility of separating 
the contributions from deformation, plastic rearrange- 
ments and large scale displacement of the block to the 
decorrelation. Nevertheless, some hints of the process 
can be inferred from a sequence before a micro-rupture, 
as shown in Figure [7J We can see on this figure a " rup- 
ture" occurring at an angle of 22.32°. Yet the failure hap- 
pens before as evidenced by a decorrelated zone mainly 
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22.14° 22.20° 22.26° 




22.32° 22.38° 22.44° 

FIG. 7. Detailed sequence describing a micro-slip event for 
an experiment with glass beads at "normal" hygrometry and 
packing fraction condition. The images shown are consecu- 
tive for angles around 22°. The angle increment between the 
states used to calculate each correlation map is 0.06°. The 
red arrows show the localisation of the failure plane. 

along a line (see arrow on fig. [7] 22.26°). This line is 
indicative of a zone of shearing between two translating 
blocks. As a matter of fact, a small solid translation does 
not modify the optical paths relative to each other and 
the images stay correlated. Nevertheless, the two blocks 
are not mechanically identical: the activity in the up- 
per section is clearly larger than in the bottom. Finally, 
the micro-rupture event occurs. It is not temporally re- 
solved and appears as a large totally decorrelated area. 
Please note that the total time of the sequence is ~ 3 s, 
which is very large compared to the inertial time corre- 
sponding to a displacement of, say, one bead diameter 
\fdfg — 5 ms. So, inertial effects do not seem to play a 
role in the development of the micro-rupture. The phe- 
nomenology stays globally the same from one experiment 
to the other, either for sand or glass beads: the density 
of rearrangements increases in the area where the rup- 
ture will take place, nevertheless a line of decorrelation 
at some depth of the sample as the one distinctly seen at 
22.26° on Figure [7] is not always clearly seen. 

We thus observe that the density of rearrangements 
and the large periodic precursors events are not indepen- 
dent. Consequently, in the following section we will focus 
on the quantitative study of the rearrangements in the 
sample. 



temporal images at different levels of the correlation, giv- 
ing an insight into the distribution of the deformation 
in the depth of the sample during the inclination. Fig- 
ure[8]shows how we proceed. The original spatiotemporal 
graph is shown on Figure Ufa). Figure E^b) and (c) are 
binarized figures obtained from Figure [8fa) for different 
values of the threshold. For Figure Efb), all values of 
correlations smaller (respectively greater) than 0.15 are 
black (resp. white). In the case of Figure [He), the pro- 
cedure is the same with a value of the threshold for the 
binarization of 0.07. We observe that during the inclina- 
tion process the level of deformation evolves as a front 
in the sample: the depth at which a fixed level of de- 
formation is reached in the sample depends linearly on 
the angle of inclination. Figure |Hf d) shows a superposi- 
tion of the fronts extracted from the binarized images for 
five different values of the thresholds. We observe that 
the slope of the front remains approximately constant. 
The solid line Figure [5(a) is obtained by averaging the 
slopes of the fronts from Figure [8f d). This shows that the 
growth of the precursors events is identical to the overall 
repartition of the deformation in the sample. Indeed, the 
average slope over several experiments, all carried out 
in "normal" conditions of packing fraction and relative 
humidity, is 13 d/° for glass beads and 8 d/° for sand. 
Those slopes are of the same order and in the same ratio 
as the values of the variation of the depth of the precur- 
sors with angle (9 d/° for glass beads and 4 d/° for sand) 
determined from data plotted on Figure [6] 




FIG. 8. (Color online) (a) Spatio-temporal representation of 
the x-averaged correlation Gi as a function of the depth and 
of the inclination angle. The material is sand in "normal" 
conditions. The dashed line gives the mean slope of the fronts 
extracted from the binarized images (see main text), (b) and 
(c) Same data binarized with two different thresholds: (b) 
0.15 and (c) 0.07. (d) Fronts extracted from the binarized 
spatio-temporal graphic for different values of the threshold, 
the gray scale is lighter when the threshold is smaller; yellow: 
0.07, green: .11, red: .15, blue: .19, black: .23. 



C. Analysis of the rearrangements 

1. Density of activity in the material 

To evidence the level of activity in the samples dur- 
ing the inclination process, we first threshold the spatio- 



In addition to this front which shows a linear relation 
between angle and depth at a value of the deformation 
and which superimposed on the growth of the precursors 
at large deformations, a lot of experiments display an- 
other deformation front very close to the surface. Such 
a phenomenon is certainly linked to a boundary effect in 
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the vicinity of the free surface. An example of such be- 
havior is shown on Figure [9j We observe in that figure 
that apart from the main front that overlaps the precur- 
sors, a thin zone of deformation breaks the slope of the 
front near the surface. 




threhold : 12% 

FIG. 9. (a) Spatio-temporal representation of the x- averaged 
correlation Gi as a function of the depth and of the inclina- 
tion angle. The material is sand in "normal" conditions, (b) 
Binarized image obtained from (a) for a threshold of 0.12. 

To conclude, the deformation in the sample appears as 
a front in a spatio-temporal representation. That front 
appears very soon in the inclination process and superim- 
posed on the precursors peaks at large deformation. The 
linear relation that connects the angle and the depth at 
which the processes are observed seems to be of the same 
nature as the linear increase of the size of the precursors 
with the angle of inclination. 



spots N occurring between z — Az/2 and z + Az/2, with 
Az the size of the slice, can be measured for different 
depths z , as a function of the current angle 6. The result 
of that image analysis is given in figure [10] where the 
number of spots at different depths is given as a function 
of the angle of inclination. Each line corresponds to a 
different depth, from 40 d under the surface to 185 d 
deep. For a given depth the counting is stopped as a 
precursor reaching that depth occur. We observe that 
at a fixed depth the number of spots increases with the 
angle while at fixed angle it decreases with the depth. At 
all the depths, the density of spots increase strongly when 
approaching the precursor event at that depth. The final 
values of the densities are of the same order at all the 
depth. 

A further analysis of those data will be done in the dis- 
cussion part. The first conclusions from the raw data are 
that the rearrangements in the sample begin as soon as 
the inclination process begins, with a decreasing density 
in the depth of the sample. The density of spots increases 
strongly before the occurrence of a precursor event at the 
corresponding depth. 



2. Spots identification 
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FIG. 10. (Color online) Number of spots N as a function 
of the angle of inclination 9 at different depths in the sam- 
ple. The average depth (in beads diameter unity d) at which 
the measurement has been made is indicated near each curve. 
The counting has been done integrating events during 1° in- 
clination and over a slice of Az ~ 37 d in depth. 

Another way to quantify the activity in the material is 
to count the number of spots, i.e. localized deformation 
in the sample, and to look at its evolution during the tilt 
process. For this, we threshold each image in a movie, 
and then we identify those spots in each image. We then 
count the number of spots that have appeared at a certain 
depth during an angle increment of 1°. The number of 



D. Influence of parameters 

In this section, we will expose the consequences of the 
modification of the humidity and packing fraction on the 
above observations. 



1. Influence of humidity 

Figure [11] shows, for sand samples with different hu- 
midity preparations, the size of the successive precursors 
as a function of the angle of inclination. The experi- 
ments have been carried out with sand prepared in the 
same manner as previously leading to a packing fraction 
around .61 for the "normal" hygrometry conditions. The 
dispersion of the results is rather large but it can be no- 
ticed that the presence of humidity tends to increase the 
number of layers of grains implied in the precursory pro- 
cess. This is in agreement with the intuitive idea that the 
presence of humidity increase the cohesion of the system 
and consequently that more grains are involved when a 
rearrangement or a precursor occurs. 

Other trends can also be noted: precursors appears 
sooner in a humid system than in a dry one. In fact, 
precursors are better defined in the humid case, as can 
be seen by comparing the spatio-temporal diagram of an 
experiment with a dry sample (Fig. UTTa)) and a humid 
sample (Fig. HTT b)). In the dry case more rearrangements 
take place giving a blurred appearance to the diagram. 
The period of the precursors is also significantly smaller 
in the dry case. 

For humid cases, a curious pairing of precursors has 
been observed several times, an example of which is 
shown in Figure [121 An oscillation between two levels 
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FIG. 11. (Color online) (a) Depth of the successive precursors 
as a function of angle for sand samples at different humidity 
preparation. Higher values of humidity mobilize more grain 
layers and the precursors begin at smaller angles, (b) and 
(c): Spatio-temporal diagrams of sand specimens prepared at 
different humidity levels, (b) Dry case (hygrometry < 20 %); 
(c) Humid case (hygrometry > 80 %). The packing fraction 
preparation of the samples is the same. 



of activity can be observed, showing an alternation be- 
tween two different regimes after each precursors. We 
also observed in some experiments an alternatation of 
two sizes for the precursors, which explain for example 
the appearance in Figure [11] of very small values of the 
precursors near the avalanche in some humid cases: these 
values alternate with much higher values. 



all. For the other half, large micro-slips event were ob- 
served with the main difference that the mobilized area 
has a clear angle with the surface of the sample instead 
of being parallel as was the case in the previously de- 
scribed micro-slips. Such a mobilized zone can be seen 
in Figures fl3T a) to (d), where the angle of the prismatic 
mobilized region to the surface is about 30°, so that the 
surface of rupture is almost horizontal in the laboratory 
frame. The mobilized zone progresses in the system as 
the angle of inclination increases. This progression oc- 
curs by regular steps at a similar period to the precursors 
obtained in samples prepared in "normal" conditions. 
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(a) 30.48° 



FIG. 13. Maps of correlation at different angles for a sample of 
sand prepared in a compacted state at "normal" hygrometry 
conditions (32 % relative humidity): (a) 30.48°, (b) 32.64° 
and (a) 32.16°. 
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FIG. 12. Spatio-temporal diagram of sand specimen prepared 
at high humidity level (above 70 %): the precursors tends to 
appear by pairs showing a different rearrangement behavior 
periodically. 



To conclude this section concerning the influence of 
preparation over the pre-avalanche phenomenology, we 
observe that the global picture of the process, i.e. the lo- 
cal rearrangements in the whole depth of the system and 
periodic micro-ruptures appearing every few degrees from 
an inclination of 15°, remains unchanged over a wide 
range of conditions of preparation. The phenomenon is 
very robust to changes in packing fraction and humidity 
conditions. Only the most compacted systems behave 
differently, presenting an internal rupture that progresses 
in the bulk of the system at an angle to the surface of 
the material. 



2. Influence of packing fraction 

When exploring different types of packing fraction 
preparation, we observed that in the case of samples pre- 
pared in the loosest state, the phenomenology is the same 
as in the "normal" conditions. In that case the activity 
in the system is very high through the depth of the pile 
and precursors are still observed. On the contrary, the 
phenomenology is modified in the most compacted piles. 
In that case, because of the significant increase of the 
avalanche angle for dense materials, we were not able to 
reach the avalanche with our experimental setup in most 
of the experiments. 

In the compacted experiments, local rearrangements 
during the inclination of the system are still observed but 
the activity in the system is much lower than in the less 
dense samples. On the contrary, large micro-failures are 
not always observed and when they are, they differ from 
the precursors events described previously. For about 
half of our experiments, no large events were observed at 



IV. DISCUSSION 

In the preceding section, we have described the com- 
plex dynamic behavior preceding the macroscopic failure 
of a granular material. The dynamics begin very early in 
the tilting process with isolated rearrangements events. 
Such events have already been reported, but with our 
sensitive side view characterization, we are able to ob- 
tain information about their spatial distribution. At the 
very early stage of the tilting process the activity is essen- 
tially limited to very superficial layers, and as the tilting 
process progresses, the activity progresses into the bulk. 
A striking feature is the presence of many micro-failure 
events in the bulk. Such events occur quite regularly 
during the tilting process. We first discuss this effect in 
the framework of Mohr- Coulomb failure criterium. After 
this, we will discuss how our observations may be related 
to the plasticity of granular material. 
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A. Cohesion 

We first discuss the localization of failure in a tilted 
granular material as it may be deduced from Coulomb 
failure criteria. The equation of static equilibrium for a 
granular material is: 



dcr x 



dx 
dcr zx 

dx 



dz 

dcTzz 

dz 



— pg sin 

— pg cos 



(2) 
(3) 



where a xz = a zx is the shear stress, a xx (a zz ) is the hor- 
izontal (vertical) normal stress, p the density, and g the 
gravity. The orientations of the axes are defined in Fig- 
ure [3j and we neglect any shear stresses in the transverse 
y direction. Assuming that the stresses are uniform in 
the x direction, we obtain: 



&zz = pg* COS t 



pgzsmt 



(4) 



The Coulomb criterion postulates that a granular ma- 
terial is stable with respect to failure if [34j : 



|r| < pa + c h , 



(5) 



for any plane inside the material. In Equation (j5]), r is the 
shear stress, a the normal stress, p the internal friction 
coefficient, and Ch the cohesion. Applying the Coulomb 
criteria to Equation (J2J), we find that the planes of failure 
are parallel to x axis, and the Coulomb criteria is now: 



tan 9 < p 



Ch 



pgz cos ( 



(6) 



First, in the absence of cohesion = 0, the depth 
of failure plane is not determined. With cohesion, the 
position of the failure plane is determined. In the case 
of constant cohesion, failure must occur at the bottom 
of the sample [35|, l3q |. where Ch/ pgz cos is the lowest. 
That result is in contradiction with the observation that 
the first micro-failure appears at a smaller depth than 
the following ones. The dependence of on the values 
of z for which failure occurs in a Mohr- Coulomb model 
is opposite to the dependence experimentally observed. 
Moreover, such a model gives only one criteria of rupture 
and is unable to predict successive micro-failures in the 
material. In any case, the Coulomb failure criterium is 
unable to explain the occurrence of the observed micro- 
failures into the bulk of the sample. 



B. Plastic deformation 

The Coulomb failure discussed in the preceding section 
does not take into account any plastic deformations be- 
fore rupture in the granular material. It is however well 
known that granular materials may yield before being 



subjected to failure. The critical state theory of soil [37j 
discussed how the plastic deformation occurs in a mate- 
rial depending on the initial state of the material. Start- 
ing from an initially loose granular material, application 
of a shear stress initially produces an elastic deformation. 
After this, yield begins with plastic deformations. Dur- 
ing this yielding, the material compacts slowly. When 
the shear stress exceeds a threshold (generally a fraction 
of the confining pressure), macroscopic failure occurs in 
the material. If the experiment begins with an initially 
dense granular material, the deformation is elastic un- 
til a maximum value of the stress is reached. When the 
stress exceeds this value, failure occurs. No plastic de- 
formation happens before that rupture and dilatancy oc- 
curs only after the failure. The fact that rearrangements 
and precursors are observed in our experiments only for 
loose or moderately dense granular samples (and not to 
densely prepared granular materials) seems to show that 
the observed deformations in our experiments correspond 
to yield before failure. We observed in a previous work, 
for a different experimental geometry, that the same kind 
of rearrangements may be associated to plastic dissipa- 
tion into the material 12711. 
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FIG. 14. (Color online) Number of spots N as a function of 
the difference between actual angle and failure angle c (z). 
Every curve corresponds to a different depth: black z = 40d, 
blue 75d, red llOd, green 150<i, yellow 185d. Inset: same 
curves in logarithmic scale. 

The link between the plastic deformations and the fail- 
ure of the material may be highlighted by considering 
the amount of plastic deformation as a function of the 
difference between the applied stress and the stress at 
which failure occurs. We have noticed (see sec lIII CTl and 
fig Ell that, at a given tilt angle (i.e. at a given applied 
stress), the density of the rearrangements decreases with 
the depth. This dependence of the density of the rear- 
rangements on the depth of the material most certainly 
originates from a gradient of properties, either due to a 
depth dependent confining pressure and/or a gradient of 
packing fraction during the preparation of the sample. 
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Such packing fraction and/or pressure dependence of the 
yield process is observed in many soil mechanics stud- 
ies, but at confining pressures noticeably larger than the 
pressure in our experiment. The tilt angle at which fail- 
ure occurs into the material (i.e. the failure stress) also 
increases with the depth (see sec JIIIB 2l and figEJ. From 
data of figinj we defined c (z) as the z-dependent value of 
the tilt angle at which failure occurs, we plot on Fig. [14] 
the number of events as a function of — c (z). The 
data obtained for different depths, which were scattered 
on figlTO]) then collapse on the same curve. This shows 
that those plastic events occur with a density which is 
given by the difference between the failure and applied 
stresses, or equivalent ly that failures happen after a cer- 
tain number of rearrangements have occurred. 

It has been shown in [27] that the rate at which lo- 
calised plastic events occur may be identified with the 
so-called fluidity of the material, which is the local rate 
of stress relaxation. The concept of fluidity has been in- 
troduced in order to explain the rheological properties of 
soft glassy material [38j. The use of this concept in or- 
der to explain rheological properties of granular material 
has also been recently proposed [39|, Hcj. The relation 
between isolated plastic events and failure may be qual- 
itatively understood. Indeed, since rearrangements may 
be seen as plastic reorganizations, we expect that each 
event redistributes some stress, accordingly to an elastic 
stress propagator [4lj |. This additional stress may also 
trigger other plastic events in their neighborhood. Be- 
cause of these processes, shear bands may form in the 
material (42j| . It follows naturally from this kind of sce- 
nario that failure occurs when a certain amount of ac- 
tivity in the material is reached. This is precisely what 
is observed: rearrangements precede micro-failure, what- 
ever is the depth. 

C. Periodicity of precursors 

Another striking observation is a rough periodicity in 
tilt angle of the precursor events. As we already shown 
(see sec lIIIB 2l and Fig. [7]), those precursors events are 
indeed failure in the bulk material. 



Failure 




FIG. 15. (Color online) Details of a slip event, (a) is the cor- 
relation image of a microslip event showing the failure plane, 
(b) is the following correlation image. The zone between the 
dashed red lines is decorrelated although below the failure 
plane, and may be interpreted as a layer of sheared granular 
material. 

We have plotted on Fig. [15] a detail of the failure al- 
ready described on Fig. [7] We see on Fig. [15]a the local- 



isation of the failure plane. One image latter (Fig. [T5lb) 
we observe that the decorrelated zone extends deeper into 
the material. So we have between the two lines of figlT5lb 
a layer of material which has been deformed and which is 
below the failure plane. We can interpret this layer as a 
zone of granular material which has been sheared by the 
motion of the upper part. This observation is in agree- 
ment with the literature about creep motion in granular 
flow. It is indeed well known that when a granular layer 
is sheared, the deformation extends into the bulk. The 
deformation decays exponentially into the bulk, with a 
characteristic length £ of the order of few grain diame- 
ters [431445] . The same decay in velocity is also observed 
during non stationary granular avalanches [46|. So the 
effect of a micro-failure event may be to produce a defor- 
mation which extends into the depth. This interpretation 
of this decorrelated layer is supported by the following 
quantitative analysis. Calling A the displacement of the 
upper block, we may expect [43-46] a displacement which 
decays as ~ A exp (— Sz/^), with Sz the distance below 
the failure plane (see fig. [T5|) . The deformation is then 
~ (A/£) exp (—5z/£). With A ~ £ c± d, we have a defor- 
mation of order 10 -6 (the limit of detection of our light 
scattering setup) for 5z ~ /n(10 6 ) x £ ~ 15d. This is 
in rough agreement with the separation between the two 
lines of flglT5lb which is ~ 15d. This agreement supports 
the hypothesis that a layer of material is sheared below 
the failure plane. 

The effect of this shear deformation may be to erase all 
the possible site for future micro-failure events because of 
the reconfiguration of the sample in that zone. The next 
micro-failure site should then be located slightly deeper 
into the bulk. The typical depth depth of this sheared 
layer slip being 15 d, we may expect from the data of 
section IIIIB 21 (variation of precursor depths with failure 
angle c± 9 d/deg for glass beads) that the angular pe- 
riod between successive micro-failure is of the order of 
~ 15/9 w 1.6 deg, a value which is close to the period 
that we measured (see figure [5]). 



V. CONCLUSION 

Our experiments on quasi-static tilting of granular 
materials show, in agreement with numerous previously 
published studies, that reorganizations occur before the 
avalanche takes place. The main advance of this study is 
to show that such rearrangements are organized spatially 
in a complex manner. For loose or moderately dense 
granular systems, we observe isolated reorganizations at 
low inclination, with a density decreasing slowly with the 
depth. As inclination increases, reorganizations occur 
under the form of micro-failure planes in the bulk, which 
are localized at increasing depths. The different local- 
izations of the failure planes imply that the underlying 
physics is different from the stick-slip already observed in 
other granular experiments [HI , where the slip plane re- 
mains the same. The micro-failures seem to occur when 
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a given level of accumulated plastic deformation is at- 
tained. This holds at every depth. The case of dense 
granular systems seems to obey to a different underlying 
physics. 

Those observations may be partially and qualitatively 
understood in the framework of yielding properties of 
granular material (e.g. Granta-Gravel material in (37| ) 
before failure. At densities below the critical one, the 
granular material yields progressively before the failure 
of the material. This yielding is continuous, and associ- 
ated to a compaction of the material. We may then in- 
terpret the observed plastic deformation as experimental 
manifestation of this yielding. Important deformations 
are then observed in the form of micro- failures, which, in 
this framework, form plastic flow in the material. The 
periodicity of such plastic flows may be understood if we 
consider collective granular flow: micro-failures involv- 
ing granular matter in the bulk with a typical extension 
length, and the next micro-failure may occur in the next 
un-deformed zone. This explanation is evidently very 
partial and only qualitative. We do not expect that a 
granular model such as Granta-Gravel will apply com- 
pletely or qualitatively, and it must be understood only 
as a starting point for future investigations. 

A puzzling point is that, to our knowledge, such col- 



lective and regular precursors motions have not been re- 
ported in numerical simulations of tilted granular mate- 
rial. It follows from our experimental results that those 
collective events may only be observed with moderately 
dense and large enough numerical granular packing. 

Another interesting point is that those micro-failures 
do not occur at the beginning of the tilt process. There 
is some minimum angle (i.e. some minimal shear stress) 
below which no micro- failures occur. To our knowledge, 
the existence of two finite critical shear stresses (one for 
micro-failure and one for the macroscopic failure) is not 
discussed in the literature. A complete understanding 
of the observation described in this paper appears as a 
challenging task. 
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